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ARTICLE INFO ABSTRACT 


Article history: 


Bovine coronavirus (BCoV) is found worldwide and causes respiratory infections and diarrhoea in calves 
and adult cattle. In order to investigate the molecular epidemiology of BCoV, 27 reverse transcription 
polymerase chain reaction (RT-PCR) positive samples from 25 cattle herds in different parts of Sweden 
were analysed. A 1038-nucleotide fragment was PCR amplified and directly sequenced. 
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vase epidemiology from outbreaks of respiratory disease or diarrhoea or from calves or adult cattle. Circulation of an iden- 
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tical BCoV strain during a 4-month period was demonstrated in calves in one dairy herd. In a regional 
epizootic of winter dysentery in Northern Sweden, highly similar BCoV strains were detected. In the 
Southern and Central regions, several genotypes of BCoV circulated contemporaneously, indicating that 
in these regions, which had a higher density of cattle than the Northern regions, more extensive trans- 
mission of the virus was occurring. Identical BCoV sequences supported the epidemiological data that 
inter-herd contact through purchased calves was important. Swedish BCoV strains unexpectedly showed 
a high homology with recently detected Italian strains. This study shows that molecular analysis of the 
spike (S) glycoprotein gene of BCoV can be a useful tool to support or rule out suspected transmission 
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routes. 


© 2012 Elsevier Ltd. All rights reserved. 


Introduction 


Bovine coronavirus (BCoV) belongs to the family Coronaviridae, 
order Nidovirales and has a single-stranded, positive-sense RNA 
genome 31 kb in length, including 13 open reading frames (Clark, 
1993; Saif, 2004). Variations in host range and tissue tropism of 
coronaviruses are attributed to the spike (S) glycoprotein (Galla- 
gher and Buchmeier, 2001). This protein is cleaved by an intracel- 
lular protease into two functional domains (Gallagher and 
Buchmeier, 2001). The S1 subunit is peripheral and responsible 
for virus binding to host-cell receptors (Kubo et al., 1994), induc- 
tion of neutralizing antibodies (Yoo and Deregt, 2001), and hae- 
magglutination activity (Schultze et al., 1991). The S1 sequence is 
variable, and mutations in this region have been associated with 
altered antigenicity and virus pathogenicity (Fazakerley et al., 
1992; Ballesteros et al., 1997). The sequence of the S2 subunit is 
more conserved and this subunit is responsible for cell membrane 
fusion activity (Yoo et al., 1991). The S1 gene has been exploited as 
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a target to study the molecular epidemiology of BCoV infection 
(Hasoksuz et al., 2002a; Jeong et al., 2005; Liu et al., 2006). 

BCoV infection is found worldwide and causes calf diarrhoea 
(CD), winter dysentery (WD) in adult cattle, and respiratory infec- 
tions in cattle of all ages (Clark, 1993; Paton et al., 1998; Saif, 
2004). Infection has a high morbidity but a low mortality (Alenius 
et al., 1991; Saif, 2004). Outbreaks typically occur in autumn and 
winter, probably in association with housing, but severe outbreaks 
in adult cattle have been reported in summer (Park et al., 2006; 
Decaro et al., 2008a). Economic losses can be heavy due to a 
marked reduction in milk yield (Travén et al., 2001; Saif, 2010). Calf 
diarrhoea and respiratory disease also cause significant losses 
(Svensson et al., 2003; Torsein et al., 2011), but the role of BCoV 
in these diseases has not been fully evaluated. 

Several serological studies have shown that BCoV infections are 
highly prevalent in Swedish cattle herds, causing disease in both 
calves and adults (Travén et al., 1999; Hagglund et al., 2006; Bido- 
khti et al., 2009). In a previous study we showed that all cattle sam- 
pled during a herd outbreak of BCoV were carrying identical BCoV 
sequences (Liu et al., 2006). The objective of the present study was 
to investigate the molecular epidemiology of BCoV by analysing 
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the polymorphic region of the S gene from herds with outbreaks of 
diarrhoea or respiratory disease and to determine if this method 
could be used to support or rule out suspected transmission routes 
between herds. 


Materials and methods 
Samples 


Faecal and nasal swab samples were collected from housed cattle in 23 dairy 
herds and two feedlot herds (Table 1) in different parts of Sweden over a 4 year per- 
iod. Samples were collected October 2005 to May 2006, October 2006 to May 2007, 
October 2007 to May 2008, and October 2008 to July 2009. Sampled cattle in all 
herds were showing clinical signs of BCoV infection and the presenting symptoms 
are summarised in Table 1. 

The old Swedish BCoV strain, C1-9202, which was detected in 1992 from a WD 
outbreak in Uppland and previously used to experimentally reproduce disease 
(Travén et al., 2001; Niskanen et al., 2002), was also sequenced in this study for 
phylogenetic comparison. Data on BCoV vaccination status was not available for 
the herds sampled in this study, but such vaccine is rarely used in Swedish cattle 
herds. 


RNA extraction, synthesis of CDNA and PCR 


RNA extraction with TRIzol LS reagent (Invitrogen) and cDNA synthesis with 
random priming was performed as described by Liu et al. (2006). Polymerase chain 
reaction (PCR) was performed as described in that paper with minor modifications. 
Briefly, 5 pL of CDNA sample was added to a tube containing 45 uL of PCR mixture, 
comprising 5 LL of 10 x PCR buffer, 1 uL of 10 mM dNTP mix, 5 pL of 1 mg/mL bo- 
vine serum albumin, 1.5 tL of each primer (10 1M), 5 LL of 25 mM MgCh, 1U of 
Taq DNA polymerase (AmpliTaq; Perkin-Elmer) and 24 lL of water. Two primers, 
S1AF and S1CR (Hasoksuz et al., 2002b) were used to amplify a 1550-nucleotide 
(nt) fragment using Biometra T3000 Thermocycler (LABREPCO). The thermocycling 
profile included 35 cycles of denaturation at 94 °C for 45 s, annealing at 50 °C for 
60 s, and extension at 72 °C for 200s, and a final extension at 72 °C for 7 min. 

Nested and semi-nested PCR (N- and SN-PCR) assays were developed to in- 
crease the sensitivity of the PCR method for clinical detection by amplifying 5 pL 
of the PCR products under the same reaction conditions described above, but with 
a 120s extension time, to amplify the target fragments using published primers 
(Hasoksuz et al., 2002b); S1AF/S1AR primer pair for the initial 654-nt fragment 
(fragment A) as SN-PCR, S1BF/S1BR primer pair for the second 490-nt fragment 
(fragment B) as N-PCR. 


Table 1 
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DNA sequencing and phylogenetic analysis 


The PCR products were purified and subsequently sequenced using the same 
primers as for PCR and an ABI PRISM BigDye Terminator v3.1 Cycle Sequencing 
Kit (Applied Biosystems) as previously described (Liu et al., 2006). Sequencing 
was performed at least twice in each direction in an ABI 3100 genetic analyser (Ap- 
plied Biosystems). The sequencing reads for each sample were assembled to a final 
1038-nt fragment of the S1 subunit spanning nucleotides 23641-24679 (amino 
acid residues 1-346 of the S glycoprotein) of BCoV strain Mebus (GenBank acces- 
sion number U00735). Multiple sequence alignment was performed with software 
BioEdit version 7.0.9.0 (Hall, 1999) with an engine based on the ClustalW algorithm. 

Maximum likelihood analysis was performed with software PHYML (Guindon 
and Gascuel, 2003), as previously described (Liu et al., 2009). The General Time 
Revisable (GTR) with substitution rate heterogeneity and proportion of invariable 
sites ‘GTR+1I+G’ was selected for inferring BCoV phylogeny, and 1000 non-para- 
metric bootstrapping analyses were performed to evaluate the robustness of the 
tree (Guindon and Gascuel, 2003). The dataset included the reference strains of 
BCoV: Mebus (GenBank accession number U00735); Quebec (AF220295); Ger- 
many/V270 (EF193075); alleged respiratory bovine coronavirus strains: LSU2 
AF058943) and OK-3 (AF058944); alleged enteric bovine coronavirus strains: F15 
D00731) and LY138 (AF058942). The analysis also included recently published se- 
quences of three Italian strains 179/07-11 (EU019216), 438/06-TN (EU814647), and 
339/06 (EF445634); two Korean strains: KWD19 (DQ389660) and KCD10 
DQ389641), and one Japanese strain: Kakegawa (AB354579). One Swedish strain 
Rc1N-99 (DQ121649) and two Danish strains Ac1F-03 (DQ121619) and Hc1F-01 
DQ121628) from the previous study (Liu et al., 2006) were included in the analysis, 
where 414 nt were treated as missing data. The analysis also included Human cor- 
onaviruses: Hu-OC43 (Z32769) and Hu-4408 (LO7748); Giraffe coronavirus 
EF424624) and Alpaca coronavirus (DQ915164). 


Results 
Genetic variation of the S gene 


Comparative analysis of the 1038 nt sequence showed that all 
the 26 Swedish strains obtained in 2005-2009 shared a high iden- 
tity both at the nt level (>97.7%) and at the deduced amino acid 
(aa) level (>96.7%). Compared with the reference strain Mebus, 
29-36 nt substitutions were identified, resulting in 16-20 aa 
changes that were clustered around residues 11-40, 88-115, 
143-188, and 241-257. 


Location of herds positive for BCoV, sampling time, animals sampled, symptoms and type of sample. 


Herd Type Year/month 
AC1 Dairy 2006/11 
AC2 Dairy 2007/11 
C1 Dairy 1992/02 
c2 Dairy 2007/01 
c3 Dairy 2007/11 
c4 Dairy 2007/12 
cS Dairy 2008/01 
C6 Feedlot 2009/03 
il Dairy 2007/02 
12 Dairy 2007/02 
13 Dairy 2007/03 
14 Dairy 2008/10 
M1 Dairy 2006/02 
M2? Dairy 2006/05 
N1 Dairy 2005/11 
N2? Dairy 2006/03 
N3 Dairy 2006/05 
N4 Dairy 2009/03 
P1 Dairy 2009/02 
U1 Feedlot 2009/07 
Y1 Dairy 2008/01 
Z1 Dairy 2007/11 
Z2 Dairy 2007/11 
Z3 Dairy 2007/11 
ZA Dairy 2007/11 
Z5 Dairy 2007/11 
Z6 Dairy 2007/12 


Age (symptoms ?) Sample type County 
Adult (WD) Faecal Vasterbotten 
Adult (WD) Faecal Vasterbotten 
Adult (WD) Faecal Uppland 
Adult (WD) Faecal Uppland 
Adult (WD) Faecal Uppland 
Adult (WD) Faecal Uppland 
Adult (WD) Faecal Uppland 
Calf (RD) Nasal Uppland 
Adult (WD) Faecal Gotland 
Adult (WD) Faecal Gotland 
Adult (WD) Faecal Gotland 
Adult (WD) Faecal Gotland 
Calf (CD) Faecal Skane 

Calf (CD) Faecal Skane 

Calf (CD) Faecal Halland 

Calf (CD) Faecal Halland 

Calf (CD) Faecal Halland 
Adult (WD) Faecal Halland 
Adult (RD) Nasal Vastergotland 
Calf (RD) Nasal Vastmanland 
Adult (WD) Faecal Angermanland 
Adult (WD) Faecal Jamtland 
Adult (WD) Faecal Jamtland 
Adult (WD) Faecal Jamtland 
Adult (WD) Faecal Jamtland 
Adult (WD) Faecal Jamtland 
Adult (WD) Faecal Jamtland 


* Herds M1and N1 visited on separate occasions, designated M2 and N 
> WD, winter dysentery; CD, calf diarrhoea; RD, respiratory disease. 
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More genetic variation was observed within Southern strains 
(23 nt changes) than within Northern and Central strains (9 nt 
changes). Sequences from five herds (Z1-Z5) in the Northern re- 
gion (Jamtland) sampled during the third year were identical. 
These sequences were also identical to the sequence U1-0907 ob- 
tained from a herd in the Central region (Vastmanland) that had an 
outbreak of respiratory disease in July 2009. The two sequences 
N1-0511 and N2-0603 obtained from the same herd with a 4- 
month interval showed 100% identity, while 6 nt changes were 
found between sequences M1-0602 and M2-0605 obtained with 
a 3-month interval from another herd. 

The BCoV strains obtained from respiratory and enteric disease 
did not show any consistent differences in the sequenced S1 re- 
gion. Strain C1-9202 shared a 97.0-98.2% identity at the nt level, 
96.2-97.4% at the deduced aa level, with the Swedish strains de- 
tected in 2005-2009. GenBank accession numbers for the BCoV se- 
quences obtained in this study are JN795141-JN795167. 


Z3-0711 
22-0711 
21-0711 
U1-0907 
25-0711 
24-0711 
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Phylogenetic analysis 


The Swedish BCoV strains obtained for this study fell into four 
clusters in the consensus tree based on maximum likelihood anal- 
ysis of the 1038-nt fragment of the S gene (Fig. 1). Cluster I in- 
cluded 20 strains from all over Sweden: three from Halland and 
one from Vdstergétland in the South; 2 from Gotland, an island 
in the east; 5 from Uppland and 1 from Vastmanland in the Central 
region; 6 from Jamtland and one each from Vdsterbotten and 
Angermanland in the North. Cluster II contained 4 Swedish strains: 
AC1-0611 from the Northern region Vasterbotten and M1-0602, 
M2-0605 and N4-0903 from the Southern regions Skane and Hal- 
land, and a Danish strain AcF-03 (Liu et al., 2006). This cluster also 
included the Italian strain Bubalus/Italy/179/07, a BCoV strain de- 
tected in buffalo calves with severe enteritis (Decaro et al., 2008b). 
Cluster III contained the two Swedish strains 12-0702 and 14-0801 
from Gotland and the strains Italy/438/06 and Italy/339/06 that 
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Fig. 1. A majority-rule consensus tree (A) of maximum likelihood analysis of the Swedish BCoV strains sampled across the country (B) during four in-stable seasons. Four 
clusters are labelled I, Il, III, and IV, respectively. Bootstrapping values (italic) are placed next to the major clades. Bar indicates changes per site. The tree is presented as rooted 
with LY-138, F15 and Hu-OC43. Herd C1 is not included on the map. 
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were isolated from an outbreak of enteric and respiratory disease 
in dairy cattle (Decaro et al., 2008a). Cluster IV contained C1- 
9202, the oldest Swedish strain from Uppland, together with the 
Swedish strain Rc1N-99 and the Danish strain HcF-01 that were 
collected in 1999 and 2001, respectively. The grouping of these 
three sequences was supported by a bootstrapping value of 62%. 


Discussion 


Twenty-seven samples from 25 herds located in different parts 
of Sweden were positive for BCoV and comparative sequence anal- 
ysis showed a high identity between the strains. There were no 
particular aa residues of the initial 346 aa of the S1 subunit that 
could differentiate strains obtained from cattle with respiratory 
disease from those obtained from calves or adult cattle with diar- 
rhoea. The results suggest that the same strains of BCoV cause nat- 
ural outbreaks of respiratory disease, CD and WD, which is in 
agreement with previous studies (Travén et al., 2001; Niskanen 
et al., 2002; Saif, 2010). 

The occurrence, in this study, of a BCoV outbreak in summer 
confirms previous reports from South Korea and Italy (Park et al., 
2006; Decaro et al., 2008a) that BCoV circulates between cattle 
herds all the year round. Furthermore, the identical 1038-nt se- 
quences of strain U1-0907, obtained from an outbreak of respira- 
tory disease in calves, and strains Z1-0711-Z5-0711, obtained 
from WD outbreaks in Jamtland, indicate that one BCoV strain 
may circulate in the cattle population for more than a year. Herd 
U1 is a feedlot which regularly purchases calves from Jamtland. 
Analysis of the whole S gene (4092 nt) of strains Z2-0711 and 
U1-0907 confirmed an identity of 99.6%, suggesting that sequence 
analysis of the S gene can be used to support or rule out suspected 
transmission routes of BCoV between regions. 

The finding of identical 1038-nt sequences of two strains ob- 
tained 4 months apart from young calves in a large dairy herd 
(N1, >200 cows) indicates, for the first time, that a BCoV strain 
may circulate for an extended period of time within large herds. 
On the other hand, strains from another large herd (M1) obtained 
3 months apart had several mutations, suggesting introduction of a 
new Strain into that herd. 

All strains obtained from Northern and Central Sweden during 
the third year were highly similar (nt identity >99.7%). These 
strains were obtained from a regional epizootic of WD, i.e. a series 
of outbreaks close in time and space. These regions were shown to 
have a low herd immunity to BCoV; 70% of the herds had BCoV 
seronegative primiparous cows in the Northern regions compared 
to 10% in the Southern parts in 2007 (Ohlson et al., 2010; A. Ohlson 
et al., unpublished data). This may be related to a lower density of 
herds in Northern regions, less transportation of animals, and/or a 
difference in biosecurity between the regions (Ohlson et al., 2010). 
The low herd immunity suggests that BCoV is not continuously cir- 
culating in the Northern parts of Sweden, and that virus causing re- 
gional epizootics is introduced from other regions. Such a 
transmission is also suggested by the finding of strain N3-0605 
in Southern Sweden prior to and highly identical (nt identity 
>99.8%) to the strain causing the Jamtland outbreaks of WD (Z1- 
Z5). 

The molecular pattern of clusters Il, II] and IV suggests that the 
BCoV strains circulating in Swedish herds are closely related to 
strains present in other parts of Europe. While it is easy to under- 
stand the clustering of strains from the Central and Southern parts 
of Sweden with strains from Denmark, which borders with Swe- 
den, grouping of Swedish strains with Italian strains was surpris- 
ing. As there were no records of animals purchased from Italy, it 
is unlikely that the virus was introduced to Sweden by transporta- 
tion of animals directly from Italy. There have, however, been ani- 


mals exported from Sweden to Central and Southern Europe in 
recent years. Samples sequenced from other European countries 
are needed in order to trace the origin and/or spread of Swedish 
strains and molecular analysis of future BCoV outbreaks from dif- 
ferent parts of the world will increase the understanding of BCoV 
epidemiology, e.g. direct and indirect transmission routes between 
herds, regions or even countries. 


Conclusions 


In this study we showed that molecular analysis of the S gene of 
BCoV was a useful tool for studying the epidemiology of this virus. 
The analysed strains showed a high sequence identity, regardless 
of whether they were obtained from outbreaks of respiratory dis- 
ease or diarrhoea or from calves or adults. Circulation of a BCoV 
strain during an extended period of time was demonstrated in 
calves in a large dairy herd. In Northern Sweden, a regional epizo- 
otic of highly similar BCoV strains was detected, while in the 
Southern and Central regions more extensive transmission of the 
virus was indicated. Finally, identical BCoV sequences supported 
epidemiological data on inter-herd contact through purchased 
calves. 
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